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, We briefly describe two statistical hadronization models, based respectively on the pres- 

f— ^ , ence and absence of light quark chemical equilibrium, used to analyze particle yields 

^s^J ■ in heavy ion collisions. We then try to distinguish between these models using K/tt 

f**-) fluctuations data. We find that while the non-equilibrium model provides an acceptable 

i description of fluctuations at top SPS and RHIC energies, both models considerably 

f— *^ , under-estimate fluctuations at low SPS energies. 

> ■ 

1 . ■ 1. Introduction 

The exploration of the thermal properties of strongly interacting matter, specifi- 
cally, of its equation of state, transport coefficients, degree of equilibration, phase 
structure, and the dependence of these on the energy and system size is one of the 
V~j ■ main objectives of heavy ion research. Thus, it is natural to try to characterize the 

soft observables in these collisions using statistical mechani cs techn iques. 

11 191^141 

C$ ' While such an approach has a long and illustrious history 1 1 | J |1 , the systematic 

and quantitative comparison of data to the statistical hadr onization (SH) model 
is a comparatively recent field. A consensus has developed I 5 | 6 | 7 | 8 | 9 | the sjj 
model can indeed fit most, if not all particle yields measured at experiments con- 
ducted at a wide range of energies. Measurements conducted at the GSI Schwerionen 
Synchrotron (SIS), BNL's Alternating Gradient Synchrotron (AGS),CERN's Super 
Proton Synchrotron (SPS), and BNL's Relativistic Heavy Ion Collider (RHIC) have 
successfully been analyzed using SH ansatze. 

However, the applied SH models differ regarding the chemical equilibration con- 
dition that is presumed. As a result, it has not as yet been possible to agree sta- 
tistical physics, if any, is responsible for the striking trends observed in the energy 
dependence of some observed hadronic yields. 

In the SH model there are two types of chemical equilibrium: all models assume 
relative chemical equilibrium^!, but some also assume absolute chemical equilibrium 
which implies the presence of just the right abundances of valance up, down, and 
strange quark pairs. If the system of produced hadrons is considered to be in ab- 
solute chemical equilibrium, then at highest heavy ion reaction energy one obtains 
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Fig. 1. (Color online) Dependence of freeze-out temperature T and baryo-chemical potential /ib on 
reaction energy in the Equilibrium (a), and non-equilibrium (b) freeze-out models. The direction 
of the arrow corresponds to increasing y^. The equilibrium dependence of T and fig in the panel 
(a) is not significantly altered by the introduction of the fitted phase space occupancy 7 S and/or 
the implementation of the Canonical ensemble for strangeness. The "star" in panel (b) corresponds 
to the point where the transition to the supercooled regime occurs and the phase space changes 
from chemically under-saturated (jq < 1) to chemically over-saturated (-y q > 1). This point also 
corresponds to the energy of the "kink" and the tip of the "horn" 



chemical freeze-out temperature T ~ 160— 170 McV. Values as low as T ~ 50 McV 
are reported at lowest reaction energies available. 

The energy dependence of the freeze-out temperature than follows the trend 
indicated in panel (a) of figure [T] as the collision energy increases, the freeze-out 
temperature increases and the baryonic density (here baryonic chemical potential 
[ib) decreases^. An increase of freeze-out temperature with y/s is expected on gen- 
eral grounds, since with increasing reaction energy a greater fraction of the energy 
is carried by mesons created in the collision, rather than pre-existing baryonsSl. 

Further refinements to the equilibriu m mo del, such as the introduction of the 
chemical non-equilibrium parameter 7 s and a canonical description of the 
system at small energies/system sizes do not materially alter the behaviour 
of temperature and chemical potential shown in the panel (a) of Fig. [T] 

When assuming chemical equilibrium, the variation in the freeze-out parameters 
T and fj>B with energy is remarkably smooth. However, there are non-continuous 
features in the energy dependence of hadronic observable s, such a s the "kink" in the 
multiplicity per number of participants and the "horn" I I I I in certain particle 
yield ratios (top panel of Fig. [2|). An effort was made to interpret this in terms 
of a shift from baryon to meson dominance ^ of the hadron yields. However, in 
the chemical equilibrium model even the simple observable like K + /tt + remains 
a smooth function of reaction energy, in contrast to the experimental results (top 
panel of Fig. [2|). Introduction of 7 S and deviations from the thermodynamic limit, 
while they help in bringing some of the model predictions closer to the data, has 
so far not managed to reproduce the sharpness of features such as the kink and the 
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Non-monotonic behaviour of particle yield ratios coul d indicate a novel reac- 
tion mechanism, e.g. onset of the deconfmement phase In such a situation, the 
smoothness of the chemical freeze-out temperature dependence on energy would be 
surprising, since it would imply that at all energies, from about 1 A GcV at SIS, to 
the highest RHIC values, there is no change in either the fireball evolution dynam- 
ics, nor any other imprint from the deconfined phase on the freeze-out condition, 
which, however is visible in the strangeness and entropy yield that K + and, respec- 
tively, 7r + represent. In particular, if the expanding system undergoes a fast conver- 
sion from a Quark Gluon Plasma (QGP) to hadrons, chemical non-equilibrium and 
super-cooling^ can be motivated by entropy and flavor conservation requirements. 

If one abandons the hypothesis of absolute chemical equilibrium, and fits phase 
space occupancies for both light and strange flavors, the systematic behaviour of 
T with energy becomes quite different as is shown in panel (b) of figure [TJ The 
two higher T values at right are for 20 (lowest SPS) and (most to right) 11.6 A 
GeV (highest AGS) reactions. In these two cases the source of particles is a hot 
chemically under-saturated (T ~ 170 MeV ) fireball. Such a system could be a 
conventional hadron gas fireball that had not the time to chemically equilibrate. 

Following the thick arrow in panel (b) of figure[T]we note that somewhat smaller 
temperatures are found with further increasing heavy ion reaction energies. Here it 
is possible ^o match the entropy of the emerging hadrons with that of a system of 
nearly massless partons when one considers supercooling to T ~ 140 MeV, while 
both light and strange quark phase space in the hadron stage acquire significant 
over-saturation with the phase space occupancy 7 g=U) d > 1 and at higher energy 
also 7 S > 1. A drastic change in the non-equilibrium condition occurs near 30 A 
GeV, corresponding to the dip point on right in panel (b) of the figure [1] (marked 
by an asterisk). At heavy ion reaction energy below (i.e. to right in panel (b) of 
figure [lj) of this point, hadrons have not reached chemical equilibrium, while at this 
point, as well as, at heavy ion reaction energy above (i.e. at and to left in panel (b) 
of figure [T]) , hadrons emerge from a much denser and chemically more saturated 
system, as would be expected were QGP formed at and above 30 A GeV. This is 
also the heavy ion reaction energy corresponding to the "kink" , which tracks the 
QGP's entropy density (higher w.r.t. a hadron gas), and the peak of the "horn" 02, 
which tracks the strangeness over entropy ratio (also higher w.r.t. a hadron gas). 

Distinguishing the two models described here by yields and particle ratios data 

has proven to be problematic. While fits with j q consistently achieve a higher 

[si 

statistical significance ™ the statistical significance of equilibrium models is well 
above the conventionally accepted minimum of 5% at all energies, so it can not be 
argued that the non-equilibrium model is "objectively better" at describing data. 
Since the two models are physically different, however, it is also not appropriate 
to just assume the equilibrium model holds because of it's greater "simplicity" 
(number of parameters). Rather, observables should be found where the two models 
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give different predictions. In the next section we will show that fluctuations of K/tt 
are an example of such observable. 

2. K/n Fluctuations in chemical (non)equilibrium 

As shown in fluctuations can be an invaluable test of light quark chemical 
non-equilibrium. This is because the scaled variance cr^ = ((AiV) 2 ) / (JV) scales 
in a different way from (N) w.r.t. temperature and the light quark phase space 
occupancy j q . 

An increase in T, in general, lowers <7jv because of the greater contribution 
of resonances that introduce particle correlations. An increase in j q , on the other 
hand, increases <7jv where N is a number carried by ir (such as the multiplicity or 
the electric charge), since the increased quantum (Bose-Einstein) contributions to 
pion yields are enhanced by the phase space over-saturation specific to j q . 

Fluctuations are in general more sensitive than yields to detector acceptance 
cuts. In addition, not understood effects, such as event-by-event volume fluctuations 
can give significant contributions to the observables that are difficult to describe 
within statistical models. The solution to second objection, which also resolves part 
of the first, is to use "dynamical" fluctuations of particle ratios 



Ni/m _ r- 2 _ , stat ^ m 

"dyn —y a Nl/N2 \ rJ Nl/N2> V 1 > 

where cr stat is the fluctuation in mixed events, that, having no resonance or quantum 
corrections, should be simply given by Poissonian scaling. uni/N2 contains such 
correlation terms, so 

*N1/N2 ~ {Nl) + m ** {Ni) {Ni) W 

where the correlation term arises due to resonances decaying into N\ and N 2 (and 
in general has to b e we ighted by an experimental reconstruction probability a). We 
refer the reader to ^1 f or details on how to compute the as described here. 

We have used the ratios and yields in the data-samples of^to fit T, (J-b.sHs, an d 
the reaction volume, in the model where j q is fitted (non-equilibrium) and wh ere 
7 g = I (chemical equilibrium). The results, together with SPS and RHIC 
experimental data are shown in Fig. [21 both for the equilibrium model (squares) 
and the non-equilibrium model (circles). We have also shown the two extreme cases 
of experimental acceptance, completely preserving the correlations (a — 1, filled 
symbols) and completely destroying them (a = 0, empty symbols). The effect of 
the correlations is not negligible, but does not alter the energy dependence of the 
fluctuations in both models. 

As can be seen, a^J™ at RHIC and the highest SPS energies is compatible with 
the non-equilibrium statistical model. Fluctuations at lower energies, however, can 
not be described by any set of chemical parameters that also describes yields. 
As expected from it's smooth variation of parameters, and smallness of quantum 
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Fig. 2. (Color online) K/tt fluctuations with model calculations. Top panels shows the ratio de- 
pendence on energy 



corrections, <7 d „ does not show a substantial variation in the equilibrium (7^ = 1 ) 
model. This model under-estimates fluctuations at all energies under consideration. 

It should be noted that these results, based on preliminary data, should be taken 
with a measure of caution. The acceptance region of NA49 is highly non-trivial (a 
is likely to be neither nor 1, and strongly momentum-dependant), and our model 
does not take this into account. While the use of dynamical fluctuations eliminates a 
some of the dependence of the results acceptance cuts, the acceptance dependence 
on correlations (for instance, K* — > Kir) remains, and introduces a systematic 
error in our model. For RHIC results, cr K ± j v ± , this correlation is negligible since 
lower-lying resonances do not produce like-sign particles in their decays. For SPS 
results, <7(if++i<--)/(7r++ Tr -), the correlation term is non-negligible but not dominant 
(~ K* /K ~ 10%), as can also be seen from the compatibility of top energy SPS 
and RHIC fluctuations. 

For p/n fluctuations the greater contribution of resonance correlations 
(A, A, E,S* and Ss all decay into pir, and form ~ 70% of the proton yields at 
all \/s), makes the statistical description of these fluctuations much less reliable 
without the inclusion of acceptance cuts. For this reason, we do not use this fluctu- 
ation in our analysis, beyond noting that both the equilibrium and non-equilibrium 
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models severely over-estimate it LLLi . We await the measurement of this fluctuation 
at RHIC, whose detectors have a very different acceptance region than SPS's NA49, 
to fully ascertain the role of limited acceptance in this fluctuation's measured value. 

We can only speculate about the reason for why fluctuations at low SPS energies 
are so above the statistical model estimate. Initial conditions have now been found, 
at RHIC, to be highly inhomogcncous. if the system produced in heavy ion colli- 
sions is perfectly thermalized throughout it's evolution ("a perfect fluid"), than the 
only trace these inhomogeneities leave at freeze-out is a fluctuation in the system 
volume, that is cancelled event by event when fluctuations of ratios are considered. 

In a system that is not perfectly equilibrated, however, "kinetic" fluctuations 
due to the random nature of each collision between the system's degrees of freedom 
can indeed arise. Perhaps this is the origin of the larger K/tt fluctuations observed 
at lower SPS energies, either due to a greater impact of the hadronic "corona" , or 
lack of equilibration in the whole system. 

To test this speculation, it is necessary to model K/tt fluctuations through 
kinetic models that include the fluctuation in the initial condition, "dynamical" 
fluctuations due to the finite number of dynamical processes that generate Kaons 
and 7rs (as well as the microscopic randomness of each process) and a model of 
the detector's acceptance cuts (as mentioned, highly non-trivial for NA49). Such a 
study goes beyond the scope of this work. 

In conclusion, we have given a description of the two statistical models cur- 
rently on the market. As a way to differentiate between these models, we have 
calculated the er^? as a function of energy, and compared the model predictions to 
experimental data. The results, while interesting, are not conclusive: The equilib- 
rium model under-estimates cr^J* at all energies, while the non-equilibrium model 
describes the higher energies SPS and RHIC acceptably. It however considerably 
under-estimates o-^J* at lower SPS energies. While this might be an indication that 
the system at these energies is not completely thermalized, and hence the random 
nature of each microscopic interaction enhances the final fluctuations beyond the 
statistical expectation, this can not for now be conclusively established. 

The author thanks the Von Humboldt foundation for the financial support given, 
Frankfurt University for the hospitality provided, as well as the QM2006 conference 
organizers for the financial support that allowed him to attend the conference. We 
also thank Michael Hauer, Mark Gorenstein, Marek Gazdzicki, Marcus Bleicher, 
Volker Koch, Sangyong Jeon, and Johann Rafelski for stimulating discussions. 
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